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ABSTRACT 







The NEPRF spectral baroclinic primitive equation with 
six layers was numerically integrated over time to examine 
the effects that vertical wind profiles have on the 
development of lee cyclogenesis. In addition, the model was 
run in both linear and nonlinear modes to isolate their 
effects on the tests. The objective was to simulate a cold 
front moving over a high mountain ridge, similar to the Alps 
or Rockies, by implementing a wind reversal profile to 
determine if this was conducive to lee 'cyclogenesis . It was 
found that the wind reversal profile produced favorable 
cyclonic growth, particularly when the model was in a linear 
mode. A nonlinear wind reversal test also produced positive 
results but only for a relatively short time; thereafter 
nonlinear interactions dampened cyclonic growth 
considerably. In addition, two tests were run that allowed 
the mountain to grow in a very short time to isolate 
inertial gravity wave interactions. The gravity waves did 
produce considerable oscillations in the two tests, but 
after 15 hours or so these two tests showed similar cyclonic 
growth to the previous tests. 
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I. 



INTRODUCTION 



Considerable interest has been directed to the study of 
leeside cyclogenesis since Eady's model (1949) described 
cyclone evolution based on a linear baroclinic mathematical 
model. Recent models with more realistic topography have 
simulated mid-latitude cyclones in the lee of large mountain 
ranges based on specific dynamic mechanisms promoting storm 
growth . 

A recent study by Smith (1986) examined leeside growth 
with a vertical wind profile similar to a cold front as it 
passes over a mountain range. Smith's "lee wave theory" 
utilized a 3-dimensional linear quasi-geostrophic model 
which predicted with partial success the time scale, 
position, size and strength of a leeside low. Based on this 
theory. Smith discovered that a leeside low would grow 
dramatically if a lid, which represents a tropopause and 
thus would allow for baroclinic instability, was introduced. 
He found that this low would grow at first by orographic 
forcing, then would grow exponentially by the baroclinic 
instability mentioned above. Thus, there were two stages of 
cyclone development. However, the initial rapid growth was 
not evident when the wind profile did not reverse with 
height. When Smith compared his theoretical results with 
actual Alpine lee cyclogenesis cases, he found good 
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qualitative but poor quantitative agreement. Some 
suggestions he proposed to improve on his work were the 
inclusion of factors such as nonlinearity, mesoscale and 
low-level blocking, and more accurate modeling of the width 
of an approaching baroclinic zone. 

Hayes (1985) investigated three dynamic mechanisms which 
previous studies indicated might be influential in mid- 
latitude leeside cyclogenesis. The three mechanisms were 
enhanced leeside baroclinic instability, continuous -mode 
growth and superposition. Two atmospheric models were used 
to study these mechanisms. The UCLA finite-difference model 
and a spectral model, developed by NEPRF by Dr. T. Rosmond, 
were implemented to study these dynamic mechanisms. The 
terrain resembled the Rocky Mountains and Hayes varies the 
width of a jet with positive shear (i.e, winds not reversing 
with height) , to simulate normal tropospheric flow over the 
Rockies. A wide jet which approximated a realistic 
baroclinic zone flowing over a high and long mountain ridge, 
with the superposition mechanism involved, gave the most 
impressive results. He concluded that the orographic 
forcing mechanism may be the catalyst needed when a 
baroclinic zone coincides with this natural leeside trough. 
It also appeared that a higher mountain ridge forced a 
deeper leeside trough to promote cyclogenesis. 

The objective of this paper is to determine if a 
vertical wind profile with a height reversal is a key factor 
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in leeside cyclogenesis. The work of Smith and Hayes formed 
the basis of this study. Chapter II describes the 
atmospheric model used in this study. It is one of the 
models used by Hayes in his work. It is the NEPRF nonlinear 
primitive equation spectral model (Rosmond) and is used in 
this research since it can be run in both a linear and 
nonlinear form. In Chapter III, the vertical wind profiles, 
the initial baroclinic disturbance and the topography are 
mathematically detailed. Chapter IV gives the results of 
the two control tests (flat terrain) and six mountain 
experiments. Two vertical wind profiles, one reversing and 
the other not reversing in height are each tested in both 
linear and nonlinear modes. The last two tests are run in a 
linear mode, but the mountain growth time was only three 
hours as compared to the 36-hour growth period for the first 
four mountain tests. Thus a comparison can be made for each 
respective wind profile based on linearity and initializa- 
tion time that might indicate peculiarities of the 
atmospheric model itself. The main variable analyzed is the 
maximum deviation from the longitudinal averaged vorticity 

I 

(Cmax) phase, so a comparison of growth rates and 
its position and movement can be made between all tests. 
While does not give a truly accurate growth factor, it 
is sufficient to determine if leeside development may occur. 
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II. MODEL DESCRIPTION 



The model used in this study is a baroclinic spectral 
transform model developed by Dr. T. Rosmond for NEPRF. A 
detailed description of this model is given by Lubeck, 
Rosmond and Williams (1977). A short summary of this model, 
implemented here, follows. 

The spectral model encompasses the nonlinear primitive 
equations for an adiabatic and hydrostatic atmosphere. 
Friction is also included, but moisture (i.e., latent heat) 
and its effects are not. The basic equations in sigma 
coordinates are as follows; 




( 2 . 1 ) 





(2.3) 



3t 



-V V0 + 




(2.4) 



da 



RT 

a 



(2.5) 
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where : 



c 


vorticity 


D 


divergence 


T 


temperature 


e 


potential temperature 


^ IT = 


surface pressure ~ 


V 


horizontal velocity vector 




geopotential height 


R 


gas constant 


Cp = 


specific heat at constant pressure 


f 


Coriolis parameter 


0 = 


vertical coordinate (o = P/ir)_ 


a “ 


vertical velocity (a = ^) 

dt 


g 


In TT 


p 


pk 


- k 


R/Cp 


F 


frictional force 


A 


longitude 


4) = 


latitude 


X 


sin (J) 



The model's prognostic variables are the vorticity and 
divergence of the wind (q/D), temperature (T) , and the 
natural log of terrain pressure (q = In Pg) . The vertical 
coordinate is defined as 
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where Pg is the surface pressure and p^ is the top of the 
model atmosphere. In this study Pt = 0 

Equations 2 . 1 through 2 . 5 are transformed into spherical 
coordinates and are represented spectrally in the 
horizontal, while finite differences are used in the 
vertical . 

The vertical structure of the model follows the 
development given by Arakawa and Suarez (1983) . In this 
study only six layers which are equally spaced in sigma are 
used (Figure 2.1). The prognostic variables listed above 
are represented spectrally in the horizontal and are 
staggered in sigma so that 5 , D, U, V, and T are carried at 
the mid-point of each layer and a is carried at the top and 
bottom of each layer. The vertical vorticity, a , is 
assumed to vanish at the upper and lower boundaries. 

Two versions of the model are used in this study. The 
linear version utilizes one wavelength in the east-west 
direction, while the nonlinear version employs three 
wavelengths. The wavenumbers when used in nonlinear 
spectral formulation are 8 for the linear version, and 0 , 8 
and 16 for the nonlinear mode. Each version is accomplished 
by forcing the time tendencies to be zero at all wavenumbers 
except for the wavenumbers employed. Cyclic continuity is 
assumed in the model. 
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Figure 2 . 1 Vertical Structure of the Spectral Model 
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The spectral formulation of the variables in the 



horizontal are represented by the following equation: 



C(A,x,a,t) 



j i c^(o,t)EeWe^ 

itt=-J n=lml 



n 



J 

z 

m=-J 



I 

n=Vl 



(2.7) 



where : 



C = some variable 



(Cn)* 



= (-1) 



-1 



-m 

=n 



m = zonal wavenumber 
= meridional index 



n 



n-|m| = the niamber of zeros between the poles 

(-1 < X < 1) of the associated Legendre 
function 



J = truncation limit 

— A = (l-l)/2 nondimensional zonal coordinate 

index (1 < 1 < 16) 



The nonlinear terms are computed using the transform 
method described by Haltiner and Williams (1980). The 
longitudinal direction is treated with a Fast Fourier 
Transform and the latitudinal direction uses Gaussian 
Quadrature. The number of latitudes, N, and longitudes, M, 
satisfy 



N 



> jJ + 1 



M 



3J + 1 



16 



The number of points are chosen so there will be aliasing 
from the product terms. For this study N = 38, and M = 16. 
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III. INITIAL CONDITIONS 



The initial wind profile and terrain are the same as 
described by Hayes (1985) . A brief description of these 
initial parameters will be discussed below. 

Hayes stated that the initial jet structure should be 
carefully selected since the vertical and horizontal 
structure of the mean tropospheric flow plays a significant 
role in the development of a cyclone. In addition, the 
size, scale and orientation of a mountain ridge can affect 
the frequency of cyclogenesis. Although the topography 
adopted from Hayes' study resembles the Rocky Mountains, it 
appears sufficient for this study since we are primarily 
concerned with a wind structure that would best promote 
cyclogenetic growth in the lee of the Alps. Those 
conditions include a cold front (i.e., winds reversing with 
height) traversing the Alpine region. 

A. MEAN WIND STRUCTURE 

The complete wind profile is expressed as: 

u[({) ,p(a , (J)) ] = -fia cos 

— — 0.5 

+ fia cos (j) [1+2 (Uy+Ug)/^ a cos (|)] ' , (3.1) 
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where a is the radius of the earth, ^ is the earth's 
rotation rate, is the upper level wind current and can 
be expressed as: 

Uu[(})/P(a, (|)) ] = Uq sech2[y(4)-4)o) ] (ln(Ps/P)/ln(Ps/Pmax) / 

(3.2) 

where Uq = 65 m/s, <i>o = 45 N, Ps = 1013.25 mb, Pmax = 2 00 mb 
and Y is the halfwidth of the jet. In this study Y = 16, 
which excludes barotropic instability. Ug, the mean surface 
current, is expressed as: 

Us((j)) = Uoosech2[ y(4)-4,q) ] (3.3) 

where Uqq = -20 m/s or +20 m/s depending on whether a wind 
reversal or a non wind reversal is implemented. The v 
component for the mean wind structure is set to zero at all 
levels. Figures 3.1 and 3.2 show the mean wind structure 
for the wind reversal and non-reversal profiles respectively 
in three dimensions. Thus we see that the wind speed varies 
linearly with the log of pressure while its latitudinal 
variation is that of a Bickley jet as described by Haltiner 
and Williams. 

B. DISTURBANCE 

The initial disturbance used in this study, which is 
adopted from Hayes (1985), is a barotropic disturbance that 
will permit us to study the effect of topography on a 
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Figure 3.1 3 -dimensional Structure of the Mean Wind 

Profile for the Wind-reversal Tests at t 
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U AT HOUR 0 



5 . 




Figure 3.2 Same as Figure 3.1, Except for Non 
Wind-reversal Tests 
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barocl inic wave . 



It consists of a weak wave which varies 



sinusoidally with longitude and has a maximum amplitude at 
45° N. The geostrophically balanced fields are expressed 
as: 



where T = 273 K, p = 1013.25 mb and A is the longitude. 
Hayes determined that wavenumber = 8 exhibited maximum 
growth and is the wavenumber used in this study for the 
linear version of the model. 

C. TERRAIN 

The topography implemented in this study is adopted from 
Hayes (1985) . It is given by 



fo[A*sin(nx) sin^ (2cj)) ] 



p ' = p $ ' /RT 



(3.4) 



u' 



-(1/foa) 9$V34) 



v' = l/(foa cos (|)) a$'/3A / T' = 0 





Zm(c|> / a) - \ 

( 0 



, 1A-A^1>4AA, 



(3.5) 
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where AA is the longitudinal grid spacing and Aq is the 
longitude where the mountain is centered. In this study 
Aq = 22.5°. Z*((f)) is expressed by 

, 4^ > (f) > 

f + 3A(1) ^ 4> > 

(3.6) 

/ 4>g > (|) > (l)g - 3A4> 

, elsewhere , 

where Zg is the mountain weight, Acf) is the latitudinal grid 
spacing, <|)jj is approximately 61.75 N and •J>s 31.25 N. 

Thus, the mountain is 22.5° wide and extends from about 73°N 
to 20°N (see Figure 3.2). 

The mountain is allowed to grow to its full height at 
time tQ. For this study's purposes tQ = 36 h for most test 
runs and tQ = 3 h for two special tests. 

Inertial gravity wave growth is reduced effectively 
using the longer growth time while the shorter growth time 
permits us to see their influence. The time increment 
scheme used to allow the mountain to reach its full height 






n, 



2 IT 

‘ ‘-w’ 2= 



Walker (1982) and others have found that smooth and rapid 
adjustments occur with this method. 

It is found in this study that the desired mountain 
height of 3000 m is truncated by the spectral parameters 
mentioned in Chapter II. Thus, the resulting height of the 
mountain in this study is 2213 m. However, it is felt this 
height is sufficient for our purposes since we are primarily 
concerned with wind profiles most conducive to cyclogenesis. 
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IV. RESULTS 



A. CONTROL RUNS 

Control runs with terrain heights set to zero over the 
entire sector are made for both the wind reversal and non- 
reversal profiles. The initial disturbance described in 
Section III.B is used in both control runs as well as the 
mountain experiments in the next section. The wind profile 
parameters for the jet structure of the control experiments 
are; 



Case A: Ug = -20 m/s, u^ = 65 m/s 

Case B: Ug = 20 m/s, u^ =65 m/s 

The primitive equation spectral model utilizes 
wavenumber 8 in its linear form since this wavenumber is 
most conducive for instability and growth (Hayes, 1985) . 

Two atmospheric levels are used in all test cases to 
determine the vertical difference in vorticity growth and 
phase speed. The levels are designated as L = 3 and L = 6, 
where L = 1/2 + a/ba and they roughly correspond to a mid- 
tropospheric height and a boundary layer top respectively. 

The maximum of the deviation from the longitudinal 

I 

averaged vorticity, ^max' phases are analyzed for 
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comparison with the mountain experiments in the next 
section. The perturbed vorticity component is observed to 
measure cyclone development while the phase is used to 
determine movement of the developing cyclonic vortex. These 
are used for comparison against the mountain experiments. 

Figures 4 . 1 and 4 . 2 show the perturbed vorticity 

_ t — 

component, ^ax' versus time and Table 4.1 lists the average 
growth of ^ax tQ to tQ + 36. The perturbed vorticity 
element for case B grows at a slightly larger rate than case 
A at the upper level (L = 3) . At the lower level the growth 
rates are nearly identical. The upper level growth is 
greater in the non wind-reversal case due to a larger mean 
flow at that level (i.e., more of a shear component of 
vorticity) . It appears that the wind reversal case has no 
advantage over the non wind-reversal case when there is no 
topography . 

The analysis of the phase speeds (Table 4.2) at 5j^ax 
reveals that case A retrogrades at the rate of -6.4 deg/day 
at L = 3 and -6.1 deg/day at L = 6. Case B results for the 
same respective levels are 27.0 and 26.2 deg/day. This 
would appear to be reasonable since the mean wind speed of 
the non-reversal case is considerably larger. This may also 
suggest that the possibility for cyclone development on the 
lee of mountains may require a mean phase speed close to 
stationarity and the wind reversal profile would fit that 
prerequisite . 
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Figure 4.1 Maximum Deviation from the Longitudinally 
Averaged Vorticity, C max ^s. Time at 
Tropospheric Level, L = 3, for Cases A and B 
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Ol-j 



Case A — ^ Case B — O 




Figure 4.2 Same as Figure 4.1, Except L = 6 
for Cases A and B 
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TABLE 4.1 





AVERAGE 


GROWTH OF 
L = 3 


^max all test runs 

L = 6 L = 


6/L = 


Case 


A 


.16 


.32 


2.00 


Case 


B 


.35 


.37 


1.06 


Test 


1 


2.51 


5.31 


2.16 


Test 


2 


.40 


.37 


.93 


Test 


3 


.42 


.94 


2.23 


Test 


4 


.48 


-1.69 


-3.52 


Test 


5 


1.80 


4.31 


2.39 


Test 


6 


.88 


.06 


.07 
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TABLE 4.2 



AVERAGE PHASE SPEED OF ^OR ALL TEST RUNS 

L = 3 L = 6 



Case 


A 


-6.34 


-6.10 


Case 


B 


27.01 


26.24 


Test 


1 


-4.86 


-6.72 


Test 


2 


28.08 


27.84 


Test 


3 


12.01 


-10.50 


Test 


4 


25.96 


24.11 


Test 


5 


-6.40 


-0.83 


Test 


6 


26.27 


28.12 
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Thus, 


overall 


inspection 


does 


not 


permit one 


to 


determine 


which 


wind profile 


is 


most 


conducive 


for 


cyclogenesis when 


no mountain 


is 


imposed on the 


mean 



atmospheric flow. 

B. MOUNTAIN TESTS 

The initial field conditions for the mountain 
experiments are the same as for the control runs described 
in Section IV. A. As mentioned previously in Chapter III, 

the mountain depicted is grossly similar in shape to the 
Rocky Mountains, but should be sufficient to use since the 
wind reversal phenomenon that is encountered in the Alpine 
Leeside situation is the main focus of this study. The 

mountain is allowed to grow to its full height (2213 m) by 
the 36th hour for the first four test runs to minimize the 
impact of the inertial gravity waves generated by the model. 
The height of the mountain was originally intended to be 
3000 m but a truncation error in the spectral representation 
produced the 2213 m height. However, it appears this height 
does not seriously degrade the output, so additional runs 
are not required. The last two test runs permit the 

mountain to grow to its full height in just three hours to 
isolate the gravity wave phenomenon. The test runs and 
their critical parameters are: 

Test 1, Linear, Wave # = 8, Wind reversal (Ug = -20 m/s, 

Uy = 65 m/s) , Mountain growth period (MGP) = 36 h. 

Test 2, Linear, Wave # = 8, Non-reversal (Ug = 20 m/s, 

Uu = 65 m/s) , MGP = 36 h. 
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Test 3, Nonlinear, Wave # = 0, 8, 16, Wind reversal, 

MGP = 36 h. 

Test 4, Nonlinear, Wave # = 0, 8, 16, Non wind-reversal, 
MGP = 36 h. 

Test 5, Linear, Wave # = 8, Wind reversal, MGP = 3 h. 

Test 6, Linear, Wave # = 8, Non wind-reversal, MGP = 3 h. 

As in the control runs, the maximum deviation from the 
longitudinally averaged vorticity, ?max' phase 

position and speed are examined. 

1. Linear Tests 

Utilizing the linear version of the primitive 
equation spectral model, tests 1, 2, 5 and 6 are run with 

the same wavenumber (8) as the control runs. The nonlinear 
versions, tests 3 and 4, utilized wavenumbers 0, 8 and 16. 
Thus the behavior of the spectral model can be explored in 
both the linear and nonlinear modes. 

Before examining the results one can assume that the 
nonlinear solutions would limit the exponential growth 
inherent in a linear solution and the nonlinear mode would 
more closely approximate a real-world scenario. In 
addition, the linearized tests would serve as a benchmark to 
other studies employing linear theory (e.g. , Smith, 1984, 
1986) . 

Figures 4.3 and 4.4 contain the maximum deviation 
from the longitudinal averaged vorticity, Cmax' versus time 
for the linear mountain experiments as well as the results 
from the two control runs. At the upper level, L = 3, Test 
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Figure 4.3 



Sane as Figure 4.1, Except L — 6 
for Cases A and B 
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Figure 4.4 Same as Figure 4.1, Except for Cases 
A and B and Tests 1 and 2 at L = 6 
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1 exhibits the largest growth from tQ to to + 36. Its 
average growth during this period is over six times larger 
than test 2 (see Table 4.1). At the lower level test 1 
experienced a growth rate 14 times larger than test 2. 
Comparing test 1 to the control runs, where the topography 
was set to zero, the comparison was even larger whereas the 
non-reversal cases exhibited no significant change at either 
tropospheric level. Thus the mountain experiment with a 
wind reversal profile, test 1, experiences the largest 
cyclonic development particularly at the top of the boundary 
layer and during the period when leeside development is 
favorable. 

The phase position of 5^^ax tests 1 and 2 and the 
control runs are shown in Figures 4.5, 4.6, 4.7 and 4.8. A 
smooth average was performed to derive a reliable phase 
speed. Table 4.2 lists these results. Important findings 
from this table are that test 1 has a retrograding speed at 
both atmospheric levels while test 2 has a steady and large 
positive phase speed. Composite results indicate that all 
of the wind reversal experiments have a retrograding 
perturbed vorticity maximum. 

Table 4.3 lists the position of Cmax both levels 
for all tests and control runs from tQ through tQ + 36. 
^max at L = 3 in test 1 is slightly behind (upstream or to 
the west of) the maximum at L = 6. Test 2, the non-reversal 
case, also has the lower level maximum ahead (west) of the 
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Figure 4 . 5 



Test 1; Phase Position of Cmax 
Time at L = 3 
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Figure 4.6 Same as Figure 4.1^ Except 
for Test 1 and L = 6 
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Figure 4.7 Same as Figure 4.1, Except for 
Test 2 and L = 3 
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Figure 4.8 Same as Figure 4.1, Except 
for Test 2 and L = 6 
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TABLE 4.3 



LONGITUDINAL POSITION OF C^ax ^ROM to TO 
to + 36 FOR ALL TEST RUNS 







to+12 


to+24 


to+36 


Case A 


L = 3 


28.1 


22.5 


42.2 


16.9 


L = 6 


30.9 


28.1 


25.3 


22.5 


Case B 
L = 3 


28.1 


42.2 


11.3 


25.3 


L = 6 


36.6 


5.6 


19.7 


33.8 


Test 1 


L = 3 


5.6 


5.6 


2.8 


0.0 


L = 6 


8.4 


8.4 


5.6 


5.6 


Test 2 


L = 3 


5.6 


42.2 


5.6 


16.9 


L = 6 


28.1 


45.0 


8.4 


2.8 


Test 3 


L = 3 


8.4 


8.4 


11.3 


16.9 


L = 6 


8.4 


8.4 


2.8 


0.0 


Test 4 


L = 3 


5.6 


42.2 


5.6 


14.1 


L = 6 


2.8 


0.0 


8.4 


2.8 


Test 5 


L = 3 


8.4 


8.4 


8.4 


2.8 


L = 6 


2.8 


5.6 


8.4 


8.4 


Test 6 


L = 3 


2.8 


5.6 


8.4 


33.8 


L = 6 


2.8 


8.4 


0.0 


45.0 
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upper level maximum except at tQ + 36 where the lower 
maximum is lagging behind the upper maximum. However, this 
seems to be an anomaly. Overall, the movement and vertical 
structure of the perturbed vorticity maxima appear to move 
with the mean wind speed in both cases. 

2 . Nonlinear Tests 

Tests 3 and 4 evaluated the nonlinear mountain 
experiments. All initial conditions are similar to tests 1 
and 2 except for the added wavenumbers 0 and 16. Figures 
4.9 and 4.10 compare the nonlinear perturbed vorticity for 
both nonlinear cases as well as the linear cases at both 
levels. Inspection of these figures as well as Table 4.1 
shows that the nonlinear wind reversal case, test 3 has 
considerably less growth than the linear case at both sigma 
levels. A vertical comparison reveals that the lower sigma 
level, L = 6, shows stronger growth rates in the linear and 
nonlinear modes for the wind reversal situation versus the 
upper level. Thus both modes show strong low-level cyclonic 
development in the first 30 h. This is not the situation 
for the non-wind reversal cases, tests 2 and 4. The growth 
rates are relatively small at both atmospheric levels. 
Overall the linear case, test 1, experiences the largest 
growth at L = 6, the lower level. 

Further inspection of Figures 4.9 and 4.10 reveals 
that the nonlinear cases, tests 3 and 4 experienced 
considerable dampening after tQ + 30 (approx. 66 h) , while 
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Figure 4.9 Same as Figure 4.1, Except for 
Tests 1, 2, 3 and 4 at L = 3 
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Figure 4.10 Same as Figure 4.1, Except for 
Tests 1, 2, 3 and 4 at L = 6 
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the linear cases increased exponentially throughout. 
Another feature observed is the undulation pattern in all 
the non wind-reversal cases. This may be due to the forced 
mountain solution imposed on the pure baroclinic solution. 

Examination of the phase speed (Figures 4.11, 4.12, 
4.13 and 4.14) reveals that test 3 at L = 3 exhibits 
stationary movement from tQ to tQ + 9, followed by a 
positive movement through tQ + 80. After tQ + 80 there 
appears to be an oscillation about the 25° longitude but 
this is difficult to verify. However, test 3 at the lower 
level has a similar negative phase speed to test 1, the 
linear wind reversal case. The non-reversal cases, tests 2 
and 4 , show comparable positive phase, speeds at both 
tropospheric levels. 

The vertical structure of C^ax shows that test 3 , 
the wind reversal case, has the upper level maximum ahead 
(east) of the lower level maximum due to the positive phase 
speed at L = 3 and a negative phase speed at L = 6 (see 
Table 4.3). This decoupling is unique in this study and 
cannot be explained here. In the non-reversal case, test 4, 
there is at least a consistent stacking of the two levels 
where the upper level trails the lower level maximum at t; 
however, due to a slightly faster phase speed at L = 3, the 
upper-level maximum moves ahead of the lower-level maximum 
by tQ + 30. 
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Figure 4.11 Test 3; Longitudinal Phase Position 
of ^ax Time at L = 3 
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Figure 4 . 12 



Same as Figure 4.11, Except for L 



6 
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Figure 4.13 



Same as Figure 4.11, Except for L 
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Figure 4.14 Same as Figure 4.11, Except for L — 6 



48 



Overall comparison reveals that through the first 36 
h both wind reversal cases, tests 1 and 3 agree 
qualitatively except for the curious meandering phase speed 
at the upper level in test 3 (the nonlinear wind reversal 
experiment) . Another significant finding is that the growth 
is more than twice as large for the wind reversal studies at 
the lower level versus the upper level, while the opposite 
is found for the non-reversal studies. 

3 . Rapid Mountain Growth Tests 

These studies have the same initial conditions of 
the previous studies except the mountain is raised to its 
full height in just three hours. The purpose is to help 
identify the large fluctuations induced by inertial gravity 
waves in the linearized spectral model. 

Inspection of the results from Figures 4.15 and 4.16 
shows rapid and relatively large oscillations of ?max 
both atmospheric levels. These oscillations appear to 
dampen to a relatively low value by the 15th hour at L = 3 
and by the 8th hour at L = 6. Cmax exhibits the most growth 
in test 5, the wind reversal case, at the lower level. This 
appears consistent with the previous wind reversal mountain 
experiments. In the non-reversal studies the short mountain 
growth period does not affect the qualitative results which 
shows that test 6 agrees quite closely with tests 2 and 4 
where the upper-level growth is greater than at the lower 
level . 



49 



VOWTICITY AMPLITUDE 




0.000001- 



20 



40 



60 

TIME (HOURS) 



80 



100 



Figure 4.15 Same as Figure 4.1^ Except for 
Tests 5 and 6 at L = 3 
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Figure 4.16 Same as Figure 4.1, Except for 
Tests 5 and 6 at L = 6 
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The phase speeds for both wind profiles are also 
compatible with the previous studies. The only deviation 
observed is in test 5 which retrogresses considerably more 
at the lower level (L = 6) . 

A vertical analysis (Table 4.3) shows consistent 
vertical stacking from tQ (3rd hour) to tQ + 36 (39th hour) 
for the wind reversal test. The non-reversal case also 
shows a similar stack through the same period, and it begins 
to move rapidly away from the mountain in unison by tQ + 39. 

From inspection of all six mountain experiments the 
most significant results ar that the wind reversal profiles 
can be twice as conducive to low-level vorticity growth 
versus upper-level growth, and the phase speeds of the wind 
reversal investigations also appear to be more favorable, 
since they are closer to stationarity (in particular test 5) 
which would allow the time necessary for cyclone 
intensification (i.e., there would be sufficient time for 
the growing vortex to interact with the basic flow) . 
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V. CONCLUSIONS 



The objective of this research is to determine whether a 
wind profile which reverses direction in the vertical can 
lead to stronger leeside cyclogenesis, than a profile which 
does not reverse. One basis of this study is the work by 
Smith (1986) who determined that a wind backing with height 
(i.e., a vertical wind profile through a cold front) is more 
conducive to cyclone growth. Smith (1986) utilized a quasi- 
geostrophic linear model. He imposed a lid which simulated 
the tropopause to allow for baroclinic instability. He 
found that a lid at 10 km with a wind reversal profile 
produced the most significant cyclonic growth. This study 
used a primitive equation baroclinic spectral model in both 
linear and nonlinear modes. This was done to see the 
computational difference of the two modes as well as to 
compare the linear results with other studies, such as Smith 
(1986). The nonlinear mode also allows an evaluation of the 
impact of wave interactions that appear in the real 
atmosphere. 

Two control tests were performed how a wind reversal and 
non wind-reversal profile behaved with no topography. Six 
mountain experiments were then tested: The first two 
(reversal and non- reversal) were in the linear mode 
(wavenumber 8) with the mountain allowed to grow to its full 
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height in 36 h. The next two tested were run with the same 
parameters as the first set, except the model was in the 
nonlinear mode (wavenumbers 0, 8, 16) . The last two tests 

were run in the linear mode but the mountain grew to its 
full height in just three hours. The last set was used to 
find out whether the inertial gravity waves in this 
primitive equation model would seriously affect the growth 
of a leeside vortex in contrast to the first four mountain 
tests . 

It was found that the most rapid growth of a leeside 
cyclone, as measured by the maximum deviation from the 
longitudinally averaged vorticity ' (5 max)/ occurred in the 
linear wind reversal case (Test 1) . The nonlinear wind 
reversal case (Test 3) grew just as rapidly in the first 12 
hours but, as expected, the nonlinear wave interactions 
dampened growth thereafter. For the non-reversal cases 
(Tests 2 and 4) , Cmax grow as quickly. In addition, 

it was observed that the growth at the lower atmospheric 
level (L = 6) was normally twice that at the higher 
atmospheric level (L = 3) for all wind reversal cases. 

This study also investigated the movement (phase speed) 
of ^max) • ma^or finding is that nearly all wind 

reversal experiments (tests 1, 3 and 5) had a relatively 

slow retrograding Cmax/ particular the short mountain 

growth experiment (Test 5) at the lower atmospheric level.. 
In contrast, the movement of Cmax non-reversal 
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cases was positive and quite fast, i.e., ?max moving 
away from the favored position for leeside development quite 
rapidly during the first 30 h. Thus the wind reversal tests 
had more favorable phase speeds during the preferred time of 
leeside cyclogenesis. The only test which did not seem to 
behave in the same manner as the others was Test 3 , the 
nonlinear wind reversal case, which showed a large positive 
phase speed at the upper level (L = 3) and thus it appeared 
that the two atmospheric levels were decoupled. All other 
tests exhibited a coupling between the lower and upper 
atmospheric perturbed vorticity maximums since the phase 
speeds at both levels were quite similar. 

Based on the results of this study, it appears there is 
good support for favorable leeside development when a wind 
reversing with height is present. However, more studies 
need to be carried out. The factors to be explored might 
include; (1) masking the orographic forced solution to 
isolate the pure perturbation field by itself; (2) a more 
realistic scenario involving a more accurate terrain model 
of the Alps; and (3) a jet structure consistent with the 
observed conditions when leeside development occurs in that 
region, e.g., southwesterly flow aloft and northwesterly 
flow at the surface rather than the artificial wind 
structure in this study. 

This study dealt with a jet structure with specific 
vertical wind profiles and a mountain ridge (2213 m high) 
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perpendicular to the initial mean flow. It examined the 
growth of leeside cyclones based on the maximum deviation 
from the longitudinal averaged vorticity. A more reasonable 
indicator of growth might be the maximum deviation from the 
time averaged vorticity. Although the mountain appeared 
high enough and long enough to include the effects of low 
level blocking, a higher ridge greater than 3000 m might be 
more realistic. The effects of a latent heat source which 
might represent the Gulf of Genoa could be included in a 
future model to see its effect on cyclone intensification. 
All these factors could possibly affect the growth of a 
leeside low. 
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